PROCEEDINGS 

OF SCIENCE 



o 
o 

(N 

> 
O 



Oh! 



> 
On 

m 

CO 
CO 



o 



Elastic 5-wave scattering phase shifts and \V ub \ from 
lattice calculations of form factors for exclusive 
semileptonic decays 



J.M. Flynn* 

University of Southampton, UK 



E-mail: j f lynnOphys . soton .ac.uk 



J. Nieves 

Universidad de Granada, Spain 



E-mail: jmnieves@ugr . es 



Omnes dispersion relations make a connection between form factors for exclusive semileptonic 
decays and phase shifts in the corresponding elastic scattering channels. We describe two appli- 
cations. In the first, we use lattice form factor calculations to learn about phase shifts in elastic 
s-wave isospin-1 /2 channels for Kn, Bn, Dk and DK scattering. The aim of the second appli- 
cation is to make the determination of the CKM matrix element magnitude \V u i,\ from exclusive 
semileptonic B — ► % decays competitive with that from inclusive decays. Here we use many sub- 
tractions in an Omnes dispersion relation to motivate a simple fitting function, allowing data to 
constrain the q 2 shape of the differential decay rate and theory, primarily lattice results, to provide 
normalisation via form factor values. 
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1. Omnes dispersion relations 

Mandelstam's hypothesis of maximum analyticity and Watson's Theorem relate the phases of 
the form factors / in exclusive semileptonic M — > % decay (where M G {B,D,K}) to the phase 
shifts in the elastic Mn —* Mn scattering amplitudes in the corresponding J p and isospin channels. 
We have 

s>s th (1.1) 



+ = T(s + ie) 2iS{s) 



f+(s-ie) T(s-ie) 

where s t h = (mM + m K ) 2 and T(s) is the scattering amplitude, related to the phase shift 8(s) by 

r(s) = i^) (e2 ' i " , - 1> <L2) 

where A is the usual kinematic function. The (inverse) scattering amplitude, in the appropriate 
isospin and angular momentum channel, is found from [1,2] 

T-\s) = -7 (s)- 1 +_L_ : J_ (1.3) 
8nciy/s th V{s) V(s th ) 

Here, V is the two-particle irreducible scattering amplitude, a is the scattering length and Iq is 
calculated from a one-loop bubble diagram. This description automatically implements elastic 



unitarity, which is necessary for the phase shift to be extracted from equation (|L2j). 
For multiple multiple subtractions, {(qf,fi) : / = 0, . . . ,«}, the Omnes result reads 



f{i) = UV <ql) *zw{h(^{ft)i\tf-m)} (14) 

7lJ Sm {s-q 1 )---(s-ql)s-q 1 

a M 2 )= n « y (^)=5, 7 , x>(4 2 )=i (i.6) 

fc=0, ^; ?y Ik k=0 

One can balance the number of subtractions against knowledge of 8. In the first application below 
we use (one or) two subtractions and form-factor input information to extract the scattering length 
in the corresponding elastic scattering channels [3]. In the second application we make many 
subtractions to motivate a simple parametrisation of the form factors for exclusive semileptonic 
B — > % decays, allowing the extraction of \V u b\ from lattice form factor results combined with 
experimental partial branching fraction information [4,5]. 

2. Elastic s-wave Kit, Biz, Dtz and DK scattering lengths 

We use lattice calculations of the scalar form factor fo(q 2 ) in exclusive semileptonic decays for 
input. In the Omnes dispersion relation we use one or two subtractions to retain dependence on the 
phase shift and apply lowest order chiral perturbation theory (ChPT) or heavy meson chiral pertur- 



bation theory (HMChPT) for the two-particle irreducible amplitudes V needed for equation (1.3). 
In our fits we can then determine the scattering length, m K a, and the form factor values, /o, at the 
chosen subtraction points. 
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Figure 1: The left hand plot shows the form factor fo(q 2 ), with a 68% error band, obtained from a fit 
using a twice-subtracted Omnes relation, implementing a linear relation between /o(0) and the scattering 
length as described in the text. Red points are the form-factor inputs and the blue square shows the result 
from [6] for /o(0) (not fitted). The right hand plot shows the isospin-1/2 Kk s-wave phase shift with a 68% 
error band (grey). The phase shift plot also shows experimental data points from [7-11]. 



2.1 Elastic s-wave Kn scattering 

For the isospin-1 /2 scalar Kn channel, the lowest order ChPT expression for V is (with f % = 
92.4 MeV) 

V(s) = ^ (m 2 K +ml - |s + ^(m| - m\ ) 2 ) . (2. 1) 

We take calculated values of the scalar form factor for decays from Nf = 2 domain wall fermion 
results by RBC [6]. Since this reference does not provide chirally-extrapopolated values for the 
form factor except at q 2 = 0, we perform our own simple chiral extrapolation, as described in [3], 
to provide input pairs (q 2 ,fo(q 2 ))- To reduce the dependence on the phase shift at large values of the 
centre-of-mass energy while retaining sensitivity to the scattering length, we use subtraction points 
at q 2 = and q\ = —0.75 GeV 2 . Our two-subtraction fit shows almost complete anticorrelation 
of /o(0) and the scattering length m n a, so we redo our fit, implementing a linear relation between 
them as a constraint (we deduce the relation from a single-subtraction fit) 1 . Our results are: 

f (q 2 ) =0.827(32), /o(0) =0.948(10), m % a = 0.179(17) (2.2) 

and our fitted form factor and phase shift are shown in Figure |]. The phase-shift plot also shows 
experimental points for comparison: we emphasise that we have not fit these data, so the agreement 
with the phase shift determined from a lattice calculation is very encouraging. Since the Omnes 
integration reaches values where massive resonance exchanges could be relevant, we estimate the 
associated uncertainties by incorporating the exchange of p and K* resonances as well as nonet 
scalar mesons with masses above lGeV, using the isopsin-1/2 Kn scattering amplitude from [12]. 
This also incorporates some next-to-leading ChPT effects. We find no appreciable changes in the 
fitted form-factor values, while the scattering length increases by 6%. We have also examined Kr] 
coupled-channel effects finding again no appreciable changes in the form-factor values and this 

'The anticorrelation is not unexpected because the lowest order ChPT expressions for/o(0) and the scattering length 
are linearly related, depending only on 1 //^ (apart from masses). 
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time a decrease of up to 5% in the scattering length. Combining these effects, we ascribe an 8% 
systematic error to the scattering length, leading to a result: 



ff%a = 0.179(17)(14). (2.3) 



2.2 Elastic s-wave Bn scattering 



For the two-particle irreducible isospin-1/2 s-wave Bn scattering amplitude we use the leading 
contact term from the heavy meson chiral perturbation theory (HMChPT) lagrangian [13], 



V{s) = ^ ( 2(mi +m£) - 3s+ ' g- V- \ . (2.4) 

We have not included a contribution from the £-channel fi*-exchange diagram depending on the 
leading HMChPT B*Bn interaction term, since this vanishes at s t h and has magnitude less than 1 % 
of that from the expression above over a large range of s. 

We take input scalar form factor values from the lattice QCD calculations by the HPQCD [14] 
and FNAL [15] collaborations, assuming that the statistical errors are uncorrected, while the sys- 
tematic errors are fully-correlated. Note that the HPQCD results are updated from those we used 
in [3], while we have also added points read off Figure 7 in [15]. We also use the lightcone sum 
rule result for / (0) = /+(0) from [16]. 

We use two subtraction points at q 2 = and g^ ax = (mg — m n ) 2 and thus perform a three- 
parameter fit to /o(0), /o(<7max) an d the scattering length m n a. We find 

/o(0) = 0.257(31), MqU = 1.18(21), m n a = 0.32(29). (2.5) 

The fitted form factor and phase shift are shown in Figure ||[ We observe that the fitted value for 
/o(^max) agrees within errors with the heavy quark effective theory prediction in the soft-pion 
limit [17], f Q (m 2 B ) = f B /f n + ^(1/mf) 1.4(2) (using f B = 189(27)MeV [18]). Our central 
phase-shift curve shows evidence for a resonance at 5.6 GeV, although we cannot give an 

upper bound for the resonance mass. 

2.3 Elastic s-wave Dn and DK scattering 



To discuss the Dn phase shift we use equation ( |2.4| ) with the obvious replacement mg — > mp. 
For the DK phase shift we project into the isospin zero channel, where the two-particle irreducible 
amplitude again takes the same form with the appropriate substitutions of masses and the replace- 
ment f K -> f K ss 1 lOMeV. 

We take input scalar form factor values from the Fermilab-MILC-HPQCD lattice QCD calcu- 
lation of reference [19]. The chiral extrapolation procedure adopted there leads to parameters for 
a Becirevic-Kaidalov (BK) [20] parametrisation of fo(q 2 ), and hence an explicit functional form, 
rather than values at a set of q 2 points. We therefore generate a toy Monte Carlo ensemble of 
BK parameters and minimise the integrated squared-difference of the BK fit-function and a twice- 
subtracted Omnes fit function to determine /o(0), /o(^max) anc ^ tne scattering length. We note that 
this fit could be avoided by using the Omnes parametrisation throughout the analysis of the lattice 
data. 
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Figure 2: B% isospin-1/2 scalar form factor and phase shift, together with 68% confidence level bounds 
(grey bands). The points on the form factor plot are the inputs from [14-16]. The dashed curves on the 
phase shift plot show the effect on the statistical uncertainty of reducing the input errors to 1 /4 of their 
current value. The intercept of the phase shift with the horizontal line at 90° indicates the position of a 
resonance. 



For the D% case, we find a scattering length m n a = 0.29(4). The output phase shift shows the 
existence of an / = 1/2 s-wave resonance at 2.2(1) GeV. 

For the DK case, we find in almost all of our Monte Carlo trials that the scattering length is 
huge, effectively infinite, telling us that ReT -1 ^) = as can be seen from equation (1.3). Hence 
there should be a resonance at threshold, (rao + m£) 2 = (2.36GeV) 2 . This can be understood 
by noting the existence of a + state, D^(2317), discovered by Babar [21], which is likely an 
isoscalar [22]. Neglecting isospin- violating decays to D^n°, this state could be considered as 
an isoscalar s-wave DK bound state. In this case, following Levinson's theorem [23], the phase 

shift close to threshold has the form n + pa-\ , where p is the centre-of-mass three-momentum. 

Three-parameter fits (two subtractions and a) show that the scattering length is effectively zero, so 
we assume that the phase shift is % over the range where the integrand of the phase-shift integral is 
significant and obtain an excellent two-parameter fit using two subtractions. 



3. \V u b\ from exclusive semileptonic B — > n decay 

For our second application we use an Omnes representation for f(q 2 ) = — q 2 )f+{q 2 ) with 
many subtractions [4, 5] to motivate the fit-function 

U{q 2 ) = -T^fl [f+(si)(ml -st)]*^ . (3.1) 

m B* 1 i=0 

We include fo information with a similar Omnes representation for f(q 2 ) = fo{q 2 ) and apply 
the constraint /+(0) = /o(0). This provides an alternative to parametrisations based on the z- 
expansion [24, 25]. Adopting the fit procedure described in [5], we combine experimental binned 
partial-branching fraction information [26-29] for / + (to determine shape) with lattice [14, 15, 
30, 31] and LCSR [16] form-factor calculations of /+ and fo (for normalisation and partial shape 
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Figure 3: Results obtained from the fit to experimental partial branching fraction data and theoretical form 
factor calculations. The top left plot shows the two form factors with their error bands, the lattice and LCSR 
input points (dots: green LCSR, red HPQCD, blue FNAL-MILC) and 'experimental' points (black triangles, 
upward-pointing for tagged and downward pointing for untagged data) constructed by plotting at the centre 
of each bin the constant form factor that would reproduce the partial branching fraction in that bin. The top 
right plot shows the differential decay rate together with the experimental inputs. The bottom plots provide 
more details of the inputs and fits by showing on the left log[(m|, — q 2 )f+ (q 2 ) as a function of q 2 , and 
on the right P(j)f+ as a function of — z [24, 25] . The dashed magenta curve in the bottom right plot is a cubic 
polynomial fit in z to the Omnes curve. 



information). From a fit with subtraction points at {0, 1/3,2/3, l}^max> we determine: 

\V ub \ = (3.47 ±0.29) x 1(T 3 
/+(0)=/o(0) = 0.245 ±0.023 
/+(?Lx/3) = 0.475 ±0.046 
/+(2^/3) = 1.07 ±0.08 

We also determine the combination \V uh \f+(0) = 8.5(8) x 10~ 4 and the total branching fraction 

B(5 ^7!T/ + v) = (1.37 ±0.08 ±0.01) x 10~ 4 (3.3) 

where the first uncertainty is from our fit and the second is from the uncertainty in the experimental 
B° lifetime. The result for \V„b\ is in striking agreement with \V u b\ extracted using all other inputs 
in CKM fits and shows some disagreement with \ V u b\ extracted from inclusive semileptonic B — > % 
decays. In Figure ^[ we show our fitted form factor and differential decay rate distribution. 



/+(<ax) =7.73 ±1.29 

= 0.338 ±0.089 
/o(2?Lx/3) = 0.520 ±0.041 
/oteLJ = 1-06 ±0.26 
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